Renewable and Sustainable Energy Reviews 32 (2014) 347-378 


Contents lists available at ScienceDirect 


Renewable and Sustainable Energy Reviews 


journal homepage: www.elsevier.com/locate/rser 


Effect of roughness geometries on heat transfer enhancement (e 
in solar thermal systems - A review 


Vipin B. Gawande *, A.S. Dhoble, D.B. Zodpe 


Department of Mechanical Engineering, Visvesvaraya National Institute of Technology, Nagpur, India 


ARTICLE INFO ABSTRACT 
Article history: The performance characteristics of a solar heater and heat exchangers can be effectively improved by 
Received 27 February 2013 using artificial roughness in different forms, shapes and sizes. Artificial roughness is provided in the form 


Received in revised form 

6 December 2013 

Accepted 4 January 2014 
Available online 1 February 2014 


of different geometries such as ribs, dimple shape roughness, wire mesh, baffles, delta winglets etc. To 
determine the effect of these geometries on thermal performance of solar heaters and heat exchangers, 
several experimental and numerical studies have been carried out by various researchers. In this paper, 
an attempt has been made to review various roughness element geometries employed in solar air heaters 
Keywords: and heat exchangers in terms of heat transfer, friction factor and flow simulation techniques. Correlations 
Solar air heater developed for heat transfer and friction factor for different roughness geometries by various investigators 


Heat exchangers in solar air heaters are presente. 
Roughness geometry 


Thermal performance 
Heat transfer 
Friction factor 


O 2014 Elsevier Ltd. AII rights reserved. 


Contents 
K aiuto ret eR) ORENSE EORR TERI TET PET "xen 348 
2... Performance analysis: of solar air liedter,.i. eocsiceisee ee p PPP Her reed CRI RT OFEN EUER EUER A m PIDE 349 
3. "Thermal performance of solar air Heater ics cens dr e or eon p Omer Deme der cx Ree I eru CREE EW YEE Oe Ca wee eee Pee 349 
34. Hydraulic performance. «ees bIEEIO IR DP YOWeeremebPpiteenpesePrsnrqeÉsbseeUbeu pex med mS Eb Ear ERES 350 
32.* -Thermo Hydraulic Performan Gey «22 arte exec mets Drew e ex POUR d OC NES SIE ven So e ea E te Fu es 350 
4... Effect of rib parameters on flow patterü. c. csisd ya to's doe ee I PIGRIG kvatrerGueremeexe*» d cu cou Ea qug Vale aoa Reba ease 350 
esses T rD"—————O»)e——————————————ÁÉÁÉÁÉÉ——— 350 
42. Effectof rib height and pitcli iier do e br erre eases Herons LOU emere p CV IO e pYu n ri eR ES VE aE 351 
43. ` Effect of inclination of rib (angle of attack) ..« «isses ore ree ere uvm ques Gan eace hh eie ri eur RR EROR RUE ec d glee bel e nre ei ae 351 
44. . Effect of V-shaping Of TiD «caosa casa yaceace ei dde a Rec eR RUE e BR E KC E ACE eR Pe E Bede d oe EC GR ce 351 
4.5. Effect of width and position of gap in continuous inclined rib... . 02... cnet eee teen eee een eeas 352 
46.  Effectof idiscretizing of V-shaped TDS: iis. occ doves sees ese ce Rianne ied chase er adea Suh Sheed RR IRR EPEE D Gud BL Eee E bonded ey deme mee aes 352 
4... Effect of riD-CEOSS SCIO soos sos eR is a Dw acta ae A levee aue Rue eene Wan Ru e or ge b ended EU Sisk Sparen Rye Sad Sew dled bee reme ede gre 352 
5. Different types of roughness geometries used in solar air heater... 2... 2... eee ne ehh hme 353 
5... -Transverse continuous TID. .-...2 oce cnc: se exe He Ra Se Rae etc eee Foe exea bled reed ele weed oa Rea ea e emus 353 
5:2.  Transverse broken TIDS-.... sce cec xxr Rec epe xe a Kex E ie DEG e NOE Hed KR ng dle Sed Fc Cede d Syed DR de qs 354 
53. Inclined continuous MDS one nece Rer Cre meRa E Rx Eee ucc xe Exe eee dee od lose reed RR ere RT SURE 354 
54. Anclined rib with gdp.iekocethRenr EC pU REC Pa DER RERUM ER Ep dU Wedoa qe d eed re ex IUS Ace REV FETU 354 
55. .Gombined inclined and transverse TID &.issici ucro pem WE et Ree quor wee ce ade Send ee e ies ead ook beau RUE 355 
ST CEDE IP EEEMT-"——————r——— 355 
5,75.. “Wedged-shaped transverse TIDS «4. verae ek CR E EAE YR e e ee D oU ever Oe e e p eto ACC RUE 357 
5.8. Combination of different rib roughness elements ............lllieeeeeeeeeee eser 357 
509. Chamfered TIDS scott enne rette Bebo PR R0 a URGE ER at Sor ab iecore alt ed p ef norte b n e wi ich EUR euadere aU Up iot. 358 
5.10. Expanded:metal mesh ribs-or wire mesbhi. 14 52s dh noh tas nte Ree ORI Rm Bobs ri e Rt ERU Cen dace id a doce ene 358 
541. Dimple/protrusion shaped geometty 5... sace re es Ree Rem UR nie at's) ce UR e ORI Eh oot CR se Je e o e RR ey bene al n ceno 359 


* Corresponding author. Tel.: +91 989 065 9696. 
E-mail address: vipingawande gmail.com (V.B. Gawande). 


1364-0321/$ - see front matter © 2014 Elsevier Ltd. All rights reserved. 
http://dx.doi.org/10.1016/j.rser.2014.01.024 


348 V.B. Gawande et al. / Renewable and Sustainable Energy Reviews 32 (2014) 347-378 


5.125 Arc Shaped TIDS e oai beso E et eren ie dett Pao Roc ORE UR UR NOR has oes SAN EA Nue ob Ee TOU ub NR x D Ra RO le eae OE DUn Salata Ca areata 360 
5:13.. Metal titi: «is ioaebooane eerte be dein Beau ace cents dub dod reto S aaNet Jedi us ose tei d ded a ng co sao M epe Re IUNII MR OK use i Rug wo 360 
5:14... W-shaped IDS «eoa oae ees prede t PRO ded OE VR tcd (o8 tese AR Ede aT lane Du EE LENS DEUR d ase ve R reae Ege bat d d pop oco 360 
5.15; Discrete W-shaped. TDS 22e ope beer dte Gob ace ctus GR ae RR Ves URINE Je uelle wire te dd C Roe cqos ape e EOD NEAR. d done tee i dee Uh 360 
5,6. Usshaped HDS; csse bos sebeeseR ese tpe Sb det Sed wi E eb ated UN a e pa REOR E e BOR SUE eae ceed cad E eR RR Gade rb qu 361 
5:17;- 3ZsShaped Tibs; usi eer Roe ee Seated te AION ec Nos £e MT IERI Eee USUS Ab abby wine SOS abu Rod e RUNE RON Mogul See abe Be 361 
518; Solid Dates PP ——ITITPMDUTIITTTESUTSTT 361 
5.19: ‘Porous and perforated bathles... srsti bee eee Over M Y e ate RI dI d erexit “Gre Oe De epe Miche WAS a Dare Vic Ure dew 366 
5:20;. Perforated Dates: PPP "--———TUTTUTMT 366 
5.21.. Delta Winslet ; «4.1 oer iR ik hea COR beth oO ORE ime I Rie uis Maeno sevens ob cece Bole sda Roe eee Atenas eA E rial 367 
5.22: Use of impinging Jet. 1er ebbe b o bei Rar dre Rr Wow C YE Y eate oS RR eee, Op RD Pur Te Magee oU ndo oe e Row 371 
6; Comipttational. analysis... saos ous once diee idR erbe cade b RR ed rer Re eor ERR e ace aoe eda ke edie rer ql d X UR REIR dee Ri dee 371 
7. Comparison of thermo hydraulic performance of roughened solar air heaters .... 0.0... cece een eere 373 
8 CONClUSIONS kense an renei a ie ee da d MEE 373 
ROT OTC OS rasp eene E a T EUER a a a a AEE nin at's E E E E E O ah dies E a A e Sa Ea e E NE NAE 377 


1. Introduction 


In the current era, when there is a continuous demand of 
energy for the economic progress and industrialization, renewable 
energy sources are playing vital role in this regard. They are used 
to design high performance heat transfer systems. Heat transfer 
enhancement in these thermal systems has numerous applications 
including cooling of electronics systems, industries, agriculture, 
space heating etc. 

Of the many alternatives, in India solar energy has high solar 
isolation. It is clean, natural and available in sufficient amount. 
Devices which are designed on the concept of utilization of this 
solar energy and works on a principle of conversion of solar 
energy into thermal energy are heat exchangers and solar air 
heaters. They are considered as the most cost effective out of all 
other designed products based on solar energy mainly in com- 
mercial and industrial applications. Schematic representation of 
the solar air heater is shown in Fig. 1. 

Solar air heaters are the devices which absorb the incoming 
solar radiations and convert it into thermal energy at the absorb- 
ing surface. The thermal efficiency of solar air heaters has been 
found to be generally poor because of low convective heat transfer 
coefficient between the absorber plate and air flowing in the duct [1]. 
This is because of the formation of laminar viscous sublayer layer 
which resist the heat transfer. The attempts adopted to enhance the 
heat transfer includes provision of artificial roughness on the under- 
side of absorber plate in the form of ribs, grooves/dimples, winglets, 
baffles, twisted tapes, mesh wires, etc. 

Artificial roughness is a passive heat transfer enhancement 
technique by which thermo hydraulic performance of a solar 
heater can be improved. As the air flows through the duct of a 
solar air heater, a laminar sub layer is formed over the observer 
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Fig. 1. Schematic diagram of conventional solar heater. 


surface that obstruct heat transfer to the flowing air, thereby 
adversely affecting the thermal performance of the solar air heater. 
In order to attain higher heat transfer coefficient, it is enviable that 
flow at the heat transferring surface is made turbulent. However 
energy for creating such turbulence has to come from fan, which 
in turn increases the power requirement. So artificially roughened 
absorber plate is considered to be a good methodology to increase 
the heat transfer coefficient since it break laminar sub layer in 
order to reduce thermal resistance. But this also causes simulta- 
neous increase in friction loss in duct. It is therefore desirable 
to create turbulence in the region very close to the heat transfer- 
ring surface i.e. in the laminar sub layer only, in order to reduce 
the friction loss with the application of artificial roughness and 
power requirement may be lessened. This can be done by keeping 
the height of the roughness elements to be small in comparison 
with the duct dimension. Fig. 2 shows the flow behavior in the 
viscous sublayer region due to presence of repeated roughness 
elements [6]. 
The objective of the present paper is 


(1) To review the research of different investigators in terms of 
roughness geometries, which can enhance convective heat 
transfer in heat exchanger and solar air heaters with minimum 
increase in friction losses [2]. 

(2) To classify various types of roughness geometries so that for a 
specific purpose particular type of geometry can be selected as 
per the need. 

(3) To discuss various types of roughness geometry in terms of 
geometrical and operating parameters. 

(4) To discuss the variation of Nusselt and friction factor depen- 
dence with respect to geometrical and operating parameters. 

(5) At the same time the various correlations developed by the 
researchers are also tabulated. The ultimate objective of the 
paper is to make available a basic platform for further research 
activities. 
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Fig. 2. Effect of roughness elements on flow field. 
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Nomenclature 

Cp specific heat of air, J/kg K 

A; surface area of absorber plate, m? 
I intensity of solar radiations, W/m? 
U, overall heat loss coefficient, W/m?K 
h heat transfer coefficient, W[m?K 
Qu useful heat gain, W 

du useful heat flux, W/m? 

m mass flow rate, kg/s 

Tpm mean plate temperature, K 

Thm mean air temperature, K 

Ti fluid inlet temperature, K 

Ta ambient temperature, K 

T5 fluid outlet temperature, K 

L length of duct, m 

li inlet length of duct, m 

L2 test length of the duct, m 

L3 outlet length of the duct, m 

k thermal conductivity of air, W/mK 
AP pressure drop, Pa 

D, Dn equivalent or hydraulic diameter of duct, m 
v velocity of air in the duct, m/s 

P pitch, m 

Pļe relative roughness pitch 

e[D relative roughness height 

glp groove position to pitch ratio 

d/w relative gap position 

gle relative gap width 

et roughness Reynolds number 

sje relative shortway length of mesh 
lje relative longway length of mesh 
I/s relative length of metal grid 

W/w relative roughness width 

e roughness height 


f friction coefficient of rough surface 
fe friction coefficient of smooth surface 
fo friction coefficient of rough surface 
BR blockage ratio 

PR baffle pitch or spacing ratio 

X stream wise Cartesian coordinate (m) 
j Colburn factor 

Re Reynolds number 

ReP* Reynolds number based on tube outside diameter 
R thermo hydraulic performance factor 
R roughness function 

g heat transfer function 


Dimensionless parameters 


Fr heat removal factor 

Nu Nusselt number 

Nu, Nusselt number for smooth surface 

Nu, Nusselt number for rough surface 

p density of air, kg/m? 

Greek symbols 

a, p angle of attack, degree 

(TA)e effective transmittance absorptance product 
o wedge angle, degree 

ô thickness of laminar sub layer, m 


Abbreviations 


VG vortex generators 
CCW counter clock wise 
CW clock wise 


2. Performance analysis of solar air heater 


Thermo hydraulic performance of a solar heater helps us to 
make an efficient system. Thermal performance related with the 
heat transfer process within the collector and hydraulic perfor- 
mance deals with the pressure drop in the duct. A conventional 
solar air heater shown in Fig. 1 is considered for detail analysis of 
thermal and hydraulic performance. Design and construction 
detail of such type of a conventional system are described by Garg 
and Prakash [3]. The simplified model as shown in Fig. 3 is consists 
of wooden rectangular box. 


* Frame - The frame is generally made of wood or sometimes 
metal. 

* Matte Black Interior - All interior surfaces are painted with a 
heat tolerant matte black paint to absorb as much of the sun’s 
heat as possible. 

* Solar absorber - This is the most vital part of the unit. The 
absorber collects the heat which is transferred to air traveling 
across and through the heated surfaces. 

* Air intake/outlet - An inlet is provided at the bottom for the 
entry of the cooler air which after picking up the heat from the 
absorber, exits the top of the unit. This happens either through a 
natural process or assisted by a thermostatically-controlled fan. 

* Glazing - Top surface of the unit is sealed with clear material 
which receives sun rays which then fall on the solar absorber 
plate and build up the interior temperature. Typical glazing 


materials are polycarbonate (Lexan or twinwall type), acrylic, or 
tempered glass. 

* Insulation - Insulation is provided on the bottom and side 
walls to avoid the heat loss from the walls. 


3. Thermal performance of solar air heater 
Hottel-Whillier-Bliss equation reported by Duffle and Beckman 


[4] is commonly used to evaluate thermal performance of solar air 
heater. It is given as 


Qu = AFa[I(za)e — Ui(Ti — Ta)] (1) 
where 

Or 
qu = Qu /As = Fo[I(za)e — UL(T; — Ta)] (2) 


The rate of useful energy gain by the flowing air through duct 


of a solar air heater may also be calculated by using the following 
solar radiations 


N 


Glazing 


Air in c C Air Out 


Absorber plate, 


Fig. 3. Simplified model of an air heater. 
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Fig. 4. Effect of rib height on laminar sub layer. 


equation: 
Qu = MCpF R(T. — Tj) = hAs(Tpm — Tim) (3) 


As discussed above, heat transfer coefficient (h) is represented 
in non-dimensional form by using relationship of Nusselt number 
(Nu) reported by Duffie and Beckman [4]. 


Ny = hL/k (4) 


Furthermore thermal efficiency of a solar air heater can be 
expressed by the following equation: 


UT; — i 


; (5) 


Nth - E I(za), — 

The above equation shows that the plot between m and 
parameter ((T; — T;)/I) can be approximated by a straight line, of 
which intercept and slope are given by the values of Fr (za). and 
Fr U; respectively. 


3.1. Hydraulic performance 


Hydraulic performance of a solar air heater concerns with 
pressure drop (AP) in the duct. Pressure drop accounts for energy 
consumption by fan to propel air through the duct. Pressure drop 
can be represented in non-dimensional form by using the following 


relationship of friction factor (f), reported by Frank and Mark [5]. 


_ (AP)D; 
— 2pLV? 


(6) 


3.2. Thermo hydraulic performance 


It is desirable that design of a collector should be made in such 
a way that it should transfer maximum heat energy to the flowing 
fluid with minimum consumption of fan energy. Therefore in 
order to analyze overall performance of a solar air heater, thermo 
hydraulic performance should be evaluated by considering ther- 
mal and hydraulic characteristics of the collector simultaneously. 
The factor is evaluated for various geometries in Section 7. 


4. Effect of rib parameters on flow pattern 
4.1. Effect of rib 


Presence of rib on the underside of the absorber plate creates 
two flow separation regions, one on each side of the rib. Separa- 
tion leads to the generation of turbulence and hence the enhance- 
ment in heat transfer as well as in the friction losses takes place. 
Fig. 4 shows the effect of rib height on laminar sub layer explained 
by Prasad et al. [7]. They reported that, 
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Fig. 5. (a) Effect of rib height and (b) effect of pitch on flow pattern. 


(1) If e « 6, roughness has no effect 


(2) If e» ó, roughness has more effect on fluid pressure as 
compared to heat transfer, due to probable interference of 
turbulence induced in the already turbulent core. 

(3) If e >ô, heat transfer has more effect with moderate effect on 
fluid pressure. 


4.2. Effect of rib height and pitch 


Flow patterns downstream of a rib with variation in rib height 
and pitch [8] are shown in Fig. 5. Reattachment of the shear layer 
does not occur for a pitch ratio of 8 because of flow separation 
downstream. In the vicinity of a reattachment point, maximum 
heat transfer occurs. Heat transfer can be similarly enhanced by 
decreasing the relative roughness pitch (P/e) for fixed relative 
roughness height (e/D) or by increasing relative roughness height 
(Ple) for fixed relative roughness pitch (e/D). Relative roughness 
pitch (P/e) can be increased up to a value of 10, beyond which 
there is a decrease in heat transfer enhancement. 


4.3. Effect of inclination of rib (angle of attack) 


Angle attack of the flow with respect to rib position i.e. 
skewness of the rib towards the flow, is another important vital 
parameter that influence the pattern of the flow. Angling of the rib 
produces span wise counter rotating secondary flows which is 
responsible for significant span wise variation of heat transfer 
coefficient. The vortices moves along the rib subsequently join the 


High Coefficient Region 


Vortices 


Low Coefficient Region 


Fig. 6. Effect of inclination of rib. 


main stream, the fluid entering near the leading end of rib and 
coming out near the trailing end as shown in Fig. 6. Turbulence 
developed by moving vortices bring in cooler channel fluid in 
contact with leading end raising heat transfer rate while the 
trailing end heat transfer is relatively lower. This effect in strong 
span wise variation of heat transfer [9]. 


4.4. Effect of V-shaping of rib 


Giving shape to a long, angled rib into v-shape helps in the 
formation of two leading ends (where heat transfer rate is high) 
and a single trailing end (where heat transfer is low) resulting in 
much large area of heat transfer. V-shaped ribs form two second- 
ary flow cells as compared to one in case of a straight angled rib 
resulting in higher overall heat transfer coefficient in case of 
v-shaped rib as shown in Fig. 7. V-shaped rib with apex facing 
downstream has a higher heat transfer as compared to that of with 
apex facing upstream [9]. 
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As the relative roughness width (W/w) is increased as a result 
of producing multiple v-ribs heat transfer further increases on 
account of formation of higher number of leading ends and 
secondary flow cells as shown in Fig. 8 [10]. According to 
researchers, there is a formation of two leading ends (where heat 
transfer rate is high) and a single trailing end (where heat transfer 
rate is low) as well as two secondary flow cells which promote 
turbulence mixing which causes increase in heat transfer. Nusselt 
number increases with increases in relative roughness width 
(W/w) and attains a maximum value corresponding to relative 
roughness width (W/w) value of 6 [11]. 


4.5. Effect of width and position of gap in continuous inclined rib 


The inclusion of gap in a rib causes secondary flow along the rib 
to join the main flow and accelerate it. This in turn helps to 
energize the retarded boundary layer flow along the surface, 
resulting in enhancement of heat transfer. Heat transfer enhance- 
ment is affected by the position of the gap with respect to leading 
and trailing edge. Position of the gap near the trailing edge, results 
in more contribution of secondary flow in energizing the main 
flow through the gap and recirculation loop in the remaining part 


Vortices 


High Coefficient Region 
- A 


Sev ~ 
C y 


High Coefficient Region 


Fig. 7. Effect of V-shaping of rib. 


Leading 

ends 
Secondary 
flow 


of the rib, thereby increasing the heat transfer rate as shown in 
Fig. 9 [12]. 

Introduction of a gap in the Multi v-shaped rib allows release of 
the secondary flow and mix with main flow through the gap as 
shown in Fig. 10 [10]. This results in its acceleration, which 
energizes the retarded boundary layer flow along the surface 
resulting in the increase of the heat transfer through the gap 
width area behind the ribs [12]. It is also reported that producing a 
gap near the leading edge does not provide sufficient strength to 
the secondary flow to energize the main flow passing through the 
gap and does not lead to increase in heat transfer. Increase in heat 
transfer is achieved by increasing the relative gap distance, 
shifting the gap towards the trailing edge. This in turn increases 
the strength of the secondary flow and thus leads in increase in 
heat transfer. 


4.6. Effect of discretizing of v-shaped ribs 


The v-shaped ribs along with staggered rib pieces in between 
further increase the number and area of heat transfer regions. 
Additional rib parameters related to the size and positioning of rib 
pieces (length ratio, B/S, segment ratio, S’/S and staggering ratio, 
P' |P) with respect to main rib produce complex interaction of 
secondary flow [13]. 


4.7. Effect of rib cross section 


The level of disturbance in the flow and the size of the 
separated region are affected by the rib cross-section. Circular 
cross section had less friction factor as compared to that of 
rectangular or square cross section ribs. This is due to reduction 
in the size of the separated region. This causes decrease in inertial 


Fig. 8. Flow pattern of secondary flow for multi v-shaped rib. 
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losses and increase in skin friction, thereby decreasing the friction 
factor. As the size of the separated region reduces, heat transfer 
gets affected critically because of decrease in the level of dis- 
turbance in flow. The reduction in the heat transfer surface area 
associated with circular cross section is the cause for Nusselt 
number decrease [14]. 


5. Different types of roughness geometries used 
in solar air heater 


Different shapes of roughness elements influencing the heat 
transfer and friction factor are discussed in the following sections. 


5.1. Transverse continuous rib 


The application of small diameter wire attached on the under- 
side of absorber plate as a artificial roughness was first introduced 
by Prasad and Mullick [15] to improve the thermal performance 
of solar air heater for drying purposes. The experimental study 
was conducted using the parameters 0.019 as relative roughness 


Top surface 


Main Flow 


RIB 


ends 


height and 12.7 as relative roughness pitch. It is investigated that 
protruding wires improve plate efficiency factor from 0.63 to 0.72 
resulting in 14% improvement in thermal performance. 

Prasad and Saini [16] investigated the effect of protrusions from 
underside of absorber surface in the form of small diameter wires 
on heat transfer and friction factor for fully developed turbulent 
flow in a solar heater duct. Experimental investigation is carried 
out using relative roughness pitch of 10, 15 and 20 and relative 
roughness height of 0.020, 0.027 and 0.033 to detect the effect of 
height and pitch of roughness on heat transfer and friction. The 
maximum value of Nusselt number and friction factor is found to 
be 2.38 and 4.25 respectively at the pitch of 10. The type and 
orientation of roughness used is shown in Fig. 11. 

Kays [17] investigated that fixing small diameter protrusion 
wires normal to flow direction on the surface of absorber plate 
helps to break laminar sub layer. It was recommended that 
protrusion wire diameter of y * =50, spaced 10-20 times diameter 
and placed within the laminar sub layer are superior to turbulence 
promoters. 

The effect of transverse wire roughness fixed on underside of 
an absorber plate is studied by Gupta et al. [18]. The experiment 
was conducted using Reynolds number in the range of 3000- 
18,000, duct aspect ratio of 6.8-11.5, relative roughness height of 
0.018-0.052 at a relative roughness pitch of 10 with a range of 
roughness Reynolds number between 5 and 70. It has been found 
that the Staton number increased initially with an increase in 
Reynolds number up to Reynolds number of 12,000 and then 
decreases for further increase in Reynolds number. 

The effect of transverse wire roughness on heat and fluid 
flow characteristics for three rectangular solar air heater ducts 
(two were roughened collectors and one was a plane surface) was 
investigated by Verma and Prasad [19]. Experiment is carried out 
using Reynolds number in the range of 5000-20,000 for high duct 
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Fig. 10. Flow pattern of secondary flow for multi v-shaped rib with gap. 
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Fig. 11. Transverse rib. 


Fig. 12. Transverse broken ribs. 


aspect ratio, relative roughness height of 0.01-0.03 at a relative 
roughness pitch of 10-40 and roughness Reynolds number range 
of 8-42. An optimum value of about 71% has been reported 
corresponding to roughness Reynolds number of 24. 


5.2. Transverse broken ribs 


Sahu and Bhagoria [20] investigated effect of 90° broken 
transverse ribs on heat and fluid flow characteristics using rough- 
ness height of 1.5 mm, duct aspect ratio value of 8, pitch in the 
range of 10-30 mm and Reynolds number in the range of 3000- 
12,000. Heat transfer enhancement was reported to be 1.25-1.4 
times over smooth duct. From the experiment maximum thermal 
efficiency was found in the order of 83.5%. Geometry used has 
been shown in Fig. 12. 


5.3. Inclined continuous ribs 


Investigators found out that inclined rib gives better heat 
transfer than transverse ribs due to generation of secondary 
flow which helps in breaking of the laminar sub layer. Han and 
Park [21] studied experimentally the effect of inclined ribs 
in narrow aspect ratio ducts which results in heat transfer 
enhancement. 

The effect of relative roughness height (e/d), inclination of rib 
with respect to flow direction and Reynolds number (Re) on the 
thermo hydraulic performance of a roughened solar air heater 
for transitionally rough flow region (5 « e* — 70) is studied by 
Gupta et al. [22]. The experimental result shows that maximum 
heat transfer and friction factor was in the order of 1.8 and 
2.7 times respectively, corresponding to values of angle of 
inclination as 60° and 70°, respectively. At a relative roughness 
(e[D) value of 0.023 and Reynolds number value of 14,000, a 
comparatively best thermo hydraulic performance was reported 
by investigator. The roughness geometry investigated has been 
shown in Fig. 13. 

An experimental investigation of forced convection heat trans- 
fer in a rectangular channel (aspect ratio AR=5) with angled 
rib turbulators, inclined at 45°, was presented by means of the 
steady-state liquid crystal thermography by Tanda [81]. To study 
the effect of rib spacing on the thermal performance of the ribbed 


Fig. 13. Inclined continuous ribs. 
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Fig. 14. Geometry of the ribbed heated surface. 


channel, the parallel ribs have been installed either onto one wall 
of the channel (1RW case) or, in-line, onto two opposite walls 
(2RW case) with a rib pitch-to-height ratio ranging from 6.66 to 
20.0. Heat transfer augmentations, relative to the smooth channel 
with the same mass flow rate, decrease with Re, as typically occurs 
in rib-roughened channels, and range from 1.6 to 2.25 for the 
1 RW case and from 1.85 to 2.55 for the 2RW case. Geometry 
investigated is shown in Fig. 14. 


5.4. Inclined rib with gap 


Aharwal et al. [23] studied the effect of width and position of 
gap in inclined split-ribs having square cross section on heat 
transfer and friction characteristics of a rectangular air heater 
duct. The increase in Nusselt number and friction factor was in 
the range of 1.48-2.59 times and 2.26-2.9 times of the smooth 
duct respectively for the range of Reynolds numbers from 3000 
to 18,000. Corresponding to a relative gap width (g/e) value of 
1.0 and relative gap position (d/W) value of 0.25, values of heat 
transfer, friction factor ratio (f/f;) and thermo hydraulic perfor- 
mance parameters were found to be maximum. It was investi- 
gated that relative gap width beyond 1.0 reduces the flow 
velocities through the gap and which reduces heat transfer as 
compared to continuous ribs. If the relative gap width was taken 
lower than 1.0, then it shrinks the passage for secondary flow 
release which reduces the turbulence intensity behind the gap 
and hence reduces heat transfer. The geometry investigated has 
been shown in Fig. 15. 

Cho et al. [24] experimentally investigated the effect of gap in 
the inclined ribs on heat transfer in square duct with rib to pitch 
height ratio of 8 and angle of attack of 60°. Experiments were 
carried out by maintaining the gap width same as the rib width 
and by varying the gap position over the duct width for parallel 
and cross rib arrangement on two opposite walls, they investi- 
gated that, the inclined rib with a downstream gap shows 
significant enhancement in heat transfer as compared to that of 
continuous inclined rib arrangement. The geometry used for study 
is shown in Fig. 16. 

An experimental study has been carried out by Kumar et al. 
[25] for enhancement of heat transfer coefficient of a solar air 
heater having roughened air duct provided with artificial rough- 
ness in the form 60° inclined discrete rib. The maximum heat 
transfer enhancement occurs for the relative roughness pitch of 
12, relative gap position of 0.35 and relative roughness height of 
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0.0498. Statistical correlations for Nusselt number and friction 
factor have been developed as a function of gap position, rib 
height (or depth), pitch and Reynolds number. These correlations 
have been found to predict the values of Nusselt number and 
friction factor with average absolute standard deviation of 3.8% 
and 3.4%, respectively. The investigated geometry has been shown 
in Fig. 17. 


5.5. Combined inclined and transverse rib 


Experimental study on heat transfer and friction characteristics 
by using a combination of inclined and transverse ribs on the 
absorber plate of solar air heater was conducted by Varun et al 
[26]. The best thermal performance was reported for relative 
roughness pitch value of 8. The investigated geometry is shown 
in Fig. 18. 


WUS 


Continuous rib 


dW =0.16 


5.6. V-shaped ribs 


Experiments were conducted by Lau et al. [27-29] to study the 
effect of discrete V-shaped ribs on turbulent heat transfer and 
friction for fully developed flow of air in a square channel. They 
found that the average Staton number for the inclined 45° and 60° 
ribs was 20-35% higher than in the 90° full rib cases. Fig. 19 shows 
the typical rib configuration. 

Han et al. [30] studied the effect of parallel and V-shaped 
staggered discrete ribs and found out that 60° staggered discrete 
V-shaped ribs provide higher heat transfer than that for the 
parallel discrete ribs. 

Momin et al. [31] studied the effect of V-shaped roughness as 
shown in Fig. 20, on heat transfer and friction characteristics of 
solar air heater duct. It was investigated that V-shape ribs with an 
angle of attack of 60° enhanced Nusselt number by 1.14 and 2.30 
times and friction factor by 2.30 and 2.83 times over inclined ribs 
and smooth plate respectively. 
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Fig. 15. Inclined ribs with gap. 
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Fig. 16. Rib arrangement with a gap. 
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The heat and fluid flow characteristics of rectangular duct 
having its one broad wall heated and roughened with periodic 
'discrete V-down rib' are experimentally investigated by Singh 
et al. [32]. As Reynolds number is increased, the Nusselt number 
increase is more for discrete V-down rib as compared to contin- 
uous V-down rib. This is probably due to increase in local Nusselt 
number downstream of gap caused by increased flow velocity 
through the gap. It is observed that at all Reynolds number, the 
Nusselt number increases with increase in relative roughness pitch 
(P/e) up to the value of 8.0 and then reduces as (P/e) is increased 
further. It is observed that at all Reynolds number, the Nusselt 
number increases as relative gap width (g/e) is increased from 
0.5 to 1.0 and thereafter it decreases as (g/e) is increased further. 
The value of Nusselt number is found to be highest at angle of 
attack (@)=60°. Also maximum value of Nusselt number is 
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Fig. 17. Geometry of 60* inclined discrete rib roughness. 
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Fig. 18. Combined inclined and transverse rib. 
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recorded at relative gap position (d/w) of 0.65. The investigated 
geometry has been shown in Fig. 21. 

Comparison of thermal performance of roughened absorber 
plate fixed with staggered discrete V-apex (up and down) was 
done by Maluwork et al. [33,34]. The increase of relative roughness 
length ratio in the range of 3-7, increases Staton number. Boosting 
in Staton number ratio was found to be of the order of 1.32-2.47. 
Rib geometry used is shown in Fig. 22. 

Karwa et al. [35] had conducted experimental study on heat 
transfer and friction in a high aspect ratio duct with transverse, 
inclined, V-up continuous and V-down continuous, V-up discrete 
ribs and V-down discrete ribs. He investigated that enhancement 
in Staton number over smooth duct was found to be 65-90%,87- 
112% 102-1372, 110-147%, 93-134%, 102-142% respectively. The 
friction factor ratio for V-up continuous and V-down continuous, 
V- up discrete ribs and V-down discrete ribs was found to be 3.92, 
3.65, 2.47 and 2.58 respectively. Roughness geometry investigated 
shown in Fig. 23. 

An experimental investigation had been carried out by Hans 
et al. [11] to study the effect of multiple v-rib roughness as shown 
in Fig. 24 on heat transfer coefficient and friction factor in an 
artificially roughened solar air heater duct. A maximum enhance- 
ment of Nusselt number and friction factor due to presence of such 
an artificial roughness has been found to be 6 and 5 times, 
respectively, in comparison to the smooth duct for the range of 
parameters considered. The maximum heat transfer enhancement 
has been found to occur for a relative roughness width (W/w) 
value of 6, while friction factor attains maximum value for relative 
roughness width (W/w) value of 10. It has also been found that 
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Fig. 20. Roughened absorber plate with V-shaped ribs. 
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Fig. 19. Typical rib configuration. 


V.B. Gawande et al. / Renewable and Sustainable Energy Reviews 32 (2014) 347-378 357 


Nusselt number and friction factor attain maxima corresponding 5.7. Wedged-shaped transverse ribs 
to angle of attack (a) value of 60° and a relative roughness pitch 


(P/e) value of 8. 


Bhagoria et al. [36] experimentally investigated the effect of 
wedge shaped transverse repeated rib roughness on one broad 


P heated wall of solar air heater duct on fluid flow characteristics. 
The experiment includes the Reynolds number range from 3000 to 


18,000, relative roughness pitch 12.12 to 60.17, relative roughness 
height 0.015-0.033. Maximum Nusselt number was obtained at a 
relative pitch of 7.57 and decreases further with the increase of 


d a pitch from 7.57 to 12.12. The friction factor keeps decreasing as the 
relative roughness pitch increases. A maximum improvement of 
i heat transfer occurs at a wedge angle of about 10° whereas Nusselt 
number decreases on either side of this wedge angle. The friction 


factor increases as the wedge angle increases. The investigated 
—— Y DF geometry has been shown in Fig. 25. 


m | 5.8. Combination of different rib roughness elements 
W 

The effect of relative roughness pitch (P/e), relative roughness 

height (e/D) and relative groove position (g/P) on the heat transfer 

coefficient and friction factor has been studied experimentally by 

Jaurker et al. [37] using rib-grooved roughness. Investigation 

— = e reveals that Nusselt number increases up to 2.7 times while 


friction factor yield up to 3.6 times as compared to smooth 


Fig. 21. Discrete V-down rib arrangement. absorber plate. The maximum enhancement in heat transfer and 


friction factor comes about for a relative roughness pitch (P/e) of 
about 6 and relative groove (g/P) value of 0.4. The higher 


turbulence intensity formed due to vortices induced in and around 
the grooves, found to be responsible for higher heat transfer rates. 
The geometry investigated has been shown in Fig. 26. 

Layek et al. [38] experimentally investigated heat transfer and 
friction factor characteristics of repeated integral transverse cham- 
Fig. 22. Discrete V-shaped ribs. fered rib groove arrangement on absorber plate of a solar air 


Transverse V — down continuous 


Flow Flow 


p 


^| 
Inclined 


V — up discrete 


V — up continuous V — down discrete 


Fig. 23. Roughness geometry. 
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Fig. 24. Roughened absorber plates. 
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Fig. 25. Wedged shaped ribs. 


heater as shown in Fig. 27. Experiments were carried out using 
four relative groove positions (g/P) values of 0.3-0.6, relative 
roughness height of 0.022-0.04, chamfer angle of 5-30* and 
Reynolds number range of 3000-21,000. The experimental result 
shows that Nusselt number and friction factor increased by 3.24 
times and 3.78 times respectively as compared to smooth duct. 


C 


5.9. Chamfered ribs 


Karwa et al. [39] experimentally studied heat transfer and 
friction characteristics for rectangular ducts, having aspect ratio 
(W/H) in the range of 4.8-12, roughened with repeated integral 
chamfered ribs as shown in Fig. 28. At chamfer angle value of 15^, 
maximum value was attained by Staton number and friction factor. 
It was studied that, Staton number decreased while friction factor 
increased as the aspect ratio (W/H) was increased. 


5.10. Expanded metal mesh ribs or wire mesh 


Saini and Saini [40] experimentally studied effect of expanded 
metal mesh as roughness geometry as shown in Fig. 29. It was 
investigated that enhancement of heat transfer coefficient and 
friction factor of the order of 4 and 5 times over the smooth duct 
corresponding to the angle of attack of 61.9° and 72° respectively 
was obtained. 

A parametric study of artificial roughness geometry of 
expanded metal mesh type in the absorber plate of solar air heater 
duct has been carried out by Gupta et al. [41] as shown in Fig. 29. 
The performance evaluation in terms of energy augmentation ratio 
(EAR), effective energy augmentation ratio (EEAR) and exergy 
augmentation ratio (EXAR) has been carried out for various values 
of Reynolds number (Re) and roughness parameters of expanded 


V.B. Gawande et al. / Renewable and Sustainable Energy Reviews 32 (2014) 347-378 359 


E— —— L—— - 


Fig. 26. (a) Rib roughness geometry and (b) rib groove geometry. 
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Fig. 27. Chamfered rib grooved artificial roughness. 
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Fig. 28. Chamfered rib geometry. 
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Fig. 29. Expanded metal mesh. 


metal mesh roughness geometry in the absorber plate of solar air 
heater duct. It was investigated that the augmentation ratios 
decrease at faster rate with Reynolds number in the order of 
EAR, EEAR and EXAR. It was also studied that augmentation ratios 
increase with increase in duct depth and intensity of solar 


Air 


Fig. 30. Dimple shaped ribs. 


Fig. 31. Dimple protrusion in staggered manner. 


radiation. It was apparent that (l/e) of 40 resulted highest EXAR 
followed by (I/e) of 55. 

Thermal performance of wire-screen metal mesh was experi- 
mentally studied by Paswan et al. [42] as artificially roughened on 
the underside of absorber plate. Thermal performance enhance- 
ment of a solar air heater was depends strongly on the diameter of 
the wire and pitch, mass flow rate, insulation and inlet tempera- 
ture. It was investigated that decreasing wire pitch, thermal 
performance of a solar air heater increases. 


5.11. Dimple/protrusion shaped geometry 


An experimental investigation was conducted by Saini et al. 
[43] using dimple shaped roughness geometry as shown in Fig. 30 
on the underside of absorber plate of the solar air heater duct. The 
maximum yield in heat transfer was found corresponding to 
relative roughness height of 0.0379 and relative roughness pitch 
of 10 and minimum value of friction factor was been found 
corresponding to relative roughness height of 0.0289 and relative 
roughness pitch of 10. 

An experimental investigation has been carried out by Bhushan 
et al. [44] for a range of system and operating parameters in order 
to analyze effect of artificial roughness on heat transfer and 
friction in solar air heater duct having protrusions as roughness 
geometry as shown in Fig. 31. Maximum enhancement of Nusselt 
number and friction factor has been found 3.8 and 2.2 times 
respectively in comparison to smooth duct for the investigated 
range of parameters. Maximum enhancement in heat transfer 
coefficient has been found to occur for relative short way length 
(S/e) of 31.25, relative long way length (L/e) of 31.25 and relative 
print diameter (d/D) of 0.294. 
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Fig. 33. Arc shape roughness. 


Fig. 34. Roughness geometry. 


An experimental investigation has been carried out by Sethi 
et al. [84] to analyze the effect of artificial roughness on heat 
transfer and friction characteristics in solar air heater duct with 
dimple shaped elements arranged in angular fashion (arc) as 
roughness elements on absorber plate as shown in Fig. 32. The 
maximum value of Nusselt number has been found corresponding 
to relative roughness height of 0.036, relative roughness pitch of 
10 and arc angle of 60°. 


5.12. Arc shaped ribs 


An experimental study has been carried out by Saini et al. [45] 
for enhancement of heat transfer coefficient of a solar air heater 
having roughened air duct provided with artificial roughness 
in the form of arc-shape parallel wire as roughness element as 
shown in Fig. 33. The maximum enhancement in Nusselt number 
has been investigated as 3.80 times corresponding to the relative 
arc angle (a@/90) of 0.3333 at relative roughness height of 0.0422. 
On the other hand, the increment in friction factor corresponding 
to these parameters has been observed 1.75 times only. 

An experimental investigation has been carried out by Yadav 
et al. [83] to study the effect of heat transfer and friction 
characteristics of turbulent flow of air passing through rectangular 
duct which is roughened by circular protrusions arranged in 
angular arc fashion as shown in Fig. 34. The investigations reported 
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Fig. 36. W-shaped rib. 


the maximum enhancement in heat transfer and friction factor as 
2.89 and 2.93 times as compared with smooth duct. 


5.13. Metal grit ribs 


Karmare et al. [46] experimentally study the effect of heat 
transfer performance of a rectangular duct with metal grit ribs as 
roughness element as shown in Fig. 35, employed on one broad 
wall, transferring heat to the air flowing through it. The effect of 
relative length, height and pitch of metal grid ribs on the heat 
transfer and friction factor has been studied for the flow range of 
Reynolds numbers 4000-17,000. It has been found that maximum 
heat transfer rate was reported for the set of parameters (l/s = 1.72, 
e[Dj — 0.044, PJe — 17.5) and maximum friction factor was observed 
for the set of parameters (/[s— 1.72, e/D,=0.044, PJe— 12.5). 
Correlation for Nusselt number and friction factor were developed 
based on the experimental results. 


5.14. W-shaped ribs 


The experimental investigation of heat transfer and friction 
factor characteristics of a rectangular duct roughened with 
W-shaped ribs arranged at an inclination with respect to the flow 
direction on its underside on one broad wall was conducted by 
Lanjewar et al. [47]. The duct had a width to height ratio (W/H) of 
8.0, relative roughness pitch (P/e) of 10, relative roughness height 
(e/Dn) of 0.03375 and angle of attack of flow (a) of 30°-75°. W ribs 
had been tested both pointing in downstream W-down and 
upstream W-up to the flow. W-down ribs provide better 
thermo-hydraulic performance than W-up and V-ribs. Highest 
enhancement in thermo-hydraulic performance for W-down ribs 
was 1.98 while it was 1.81 for W-up ribs in the range of parameters 
investigated. The geometry has been shown in Fig. 36. 


5.15. Discrete W-shaped ribs 


Kumar et al. [48] experimentally investigated heat transfer and 
friction characteristics in a solar air heater having absorber plate 
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Fig. 37. Discrete W shaped ribs. 
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Fig. 39. Z-shaped rib. 


roughened with discrete W-shaped ribs. It was investigated that, 
thermal performance of the roughened channel was 1.2-1.8 times 
the smooth channel for range of parameters investigated. The 
maximum enhancement of Nusselt number and friction factor was 
reported as 2.16 times and 2.75 times that of smooth duct, 
analogous to angle of attack of 60° and relative roughness height 
of 0.0338. The geometry investigated has been shown in Fig. 37. 


5.16. U-shaped ribs 


An experimental investigation had been carried out by Bopche 
et al. [49] to study the heat transfer coefficient and friction factor 
by using artificial roughness in the form of inverted U-shaped 
turbulators as shown in Fig. 38, on the absorber surface of an air 
heater duct. As compared to the smooth duct, the enhancement in 
heat transfer and friction factor was reported of the order of 2.82 
and 3.72 times, respectively. At Reynolds number, Re— 3800, the 
maximum enhancement in Nusselt number and friction factor are 
of the order of 2.388 and 2.50, respectively. Better thermo 
hydraulic performance was reported due to generation of turbu- 
lence only in the viscous sub layer region. 
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Fig. 40. Baffle arrangement on absorber plate. 
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Fig. 41. Staggered diamond shape baffle. 


5.17. Z-shaped ribs 


Sriromrein et al. [50] studied heat transfer and friction char- 
acteristic of a rectangular duct roughened artificially with Z-shape 
ribs as shown in Fig. 39. The enhancement in heat transfer rate 
and best thermal performance was reported for Z-rib inclined at 
45°. 


5.18. Solid baffles 


The efficiency of solar air heaters with baffles as shown in 
Fig. 40 was experimentally investigated by Yeh et al. [51]. Solar air 
heaters with fins in the collector were provided with attached 
baffles which creates air turbulence and an extended heat transfer 
area. This arrangement caused a considerable improvement in 
collector efficiency. It was studied that increasing the density of 
baffles also increases the collector efficiency but this will increase 
the power consumption by blower. 

Theoretical investigation of the effect of collector aspect ratio 
on collector efficiency of baffled solar air heaters was carried out 
by Yeh et al. [52]. They investigated that with increase in aspect 
ratio, the collector efficiency is also increases. Numerical investi- 
gation of laminar periodic fluid flow and heat transfer in the 
entrance region of a two dimensional horizontal channel with 
isothermal walls and with staggered baffles was conducted by 
Sripattanapipat et al. [53] as shown in Fig. 41. Effects of different 
baffle tip angles on heat transfer and pressure loss in the channel 
were studied and the results of the diamond baffle were compared 
with those of the flat baffle. The order of enhancement was 
reported about 200-680% for using the diamond baffles. However, 
as predicted, the augmentation was associated with enlarged 
friction loss ranging from 20 to 220 times above the smooth 
channel. 

Numerical simulations of three-dimensional laminar forced 
convection flow adjacent to backward-facing step in rectangular 
duct were presented by Nie et al. [54] to inspect effects of the 
baffle on flow and heat transfer distributions. Investigators main- 
tained step height as constant and baffle were mounted onto the 
upper wall and its distance from the backward-facing step is 
varied. They reported that this arrangement of baffles increases 
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Fig. 42. Schematic of computational domain. 


Fig. 43. Multiple 60° V-baffles. 


Fig. 44. Computational domain of periodic flow. 


the magnitude of maximum Nusselt number at the stepped wall. 
The maximum Nusselt number on the stepped wall was developed 
near the sidewall, and it moved further downstream as the location 
of the baffle moves in the streamwise direction. The friction 
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Fig. 45. Duct geometry and computational domain with baffles. 
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Fig. 46. Duct geometry and computational domain. 


coefficient at the stepped wall decreases as the distance of the 
baffle from the inlet increases. Schematic of computational domain 
is shown in Fig. 42. 

Promvonge [55] carried out an experimental investigation to 
assess turbulent forced convection heat transfer and friction loss 
behaviors for airflow through a channel fitted with a multiple 60° 
V-baffle turbulator as shown in Fig. 43. Experiment was carried 
out using a channel of aspect ratio, AR= 10 and height, H=30 mm 
with three different baffle blockage ratios, (e/H=0.10, 0.20 and 
0.30) and three baffle pitch spacing ratios, (PR=P/H=1, 2 and 3) 
while the transverse pitch of the V-baffle is set to 2H and kept 
constant. 

The experimental results confirmed that the V-baffle provides 
the drastic increase in Nusselt number, friction factor and thermal 
enhancement factor values over the smooth wall channel due to 
better flow mixing from the formation of secondary flows induced 
by vortex flows generated by the V-baffle. Significant increase in 
Nusselt number and friction factor values were found for the rise 
in blockage ratio and/or for the decrease in pitch ratio values. At 
lower Reynolds number, the use of V-baffle with e/H=0.10 and 
PR=1, reported the highest value of 1.87 for thermal enhancement 
factor. 
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Fig. 47. Computational domain with baffles. 


Promvonge et al. [56] carried out a numerical investigation on 
periodic laminar flow and heat transfer behaviors in a three- 
dimensional isothermal wall square-channel fitted with 30° 
angled baffles on two opposite channel walls as shown in 
Fig. 44. The angled baffles with the attack angle of 30° were 
mounted periodically and inline arrangement on the lower and 
upper channel walls. This arrangement caused a pair of stream 
wise counter-rotating vortex (P-vortex) through the channel. They 
studied the effects of different baffle heights and three pitch ratios 
on heat transfer and flow behaviors in the channel. 

Heat transfer enhancement was reported as 1.2-11 times for 
using the 30° baffle pair with BR=0.1-0.3 and PR=1, 1.5 and 2. 
Penalty of pressure loss in the range 2 to 54 times above the 
smooth channel was also reported by investigators. The maximum 
value of thermal enhancement factor was found to be 4.0 at 
highest Reynolds number and PR=2. When effect of the baffle PR 
and BR values on heat transfer rate is examined, it was found that 


the maximum enhancement factor is about 4.0 for BR=0.15, 
PR=2.0 and Re=2000. 

Promvonge et al [57] carried out numerical study on periodic 
laminar flow and heat transfer behaviors in a three dimensional 
isothermal wall square channel fitted with 30° angled baffles on 
one channel wall as shown in Fig. 45. The range of parameters 
used for the experiment were Reynolds number in the range of 
100-2000, PR=1,1.5 and 2, BR=0.1-0.5. In this study, heat transfer 
enhancement was reported as 1.00-9.23 times for using baffle 
with BR=0.1-0.5. The pressure loss penalty was reported in the 
range of 1.09 to 45.31 times over the smooth duct. Maximum value 
of thermal enhancement factor was found to be 3.1 at BR- 0.3 and 
PR— 1.5. 

A numerical investigation had been carried out by Promvonge 
et al. [58] to examine periodic laminar flow and heat transfer 
characteristics in a three-dimensional isothermal wall channel 
of aspect ratio, AR=2 with 45° staggered V-baffles as shown in 
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Fig. 46. V-baffles with an attack angle of 45^ were mounted in 
tandem and staggered arrangement on the lower and upper walls 
of the channel to generate two pair of main stream wise vortex 
flows through the tested section. Heat transfer enhancement was 
about 100-1100% for using both the 45° baffles with BR=0.05-0.3 
and pressure loss ranging from 2 to 90 times above the smooth 
channel. Maximum value of thermal enhancement factor was 
reported as about 2.6-2.75 indicating higher thermal performance 
over smooth channel. 

A numerical investigation on laminar flow and heat transfer 
characteristics in a three dimensional isothermal wall square- 
channel fitted with inline 45^ V-shaped baffles on two opposite 
wall was done by Promvonge et al. [59]. The inline V-baffles with 
its V-tip pointing downstream and the attack angle (or half V-apex 
angle) of 45° relative to the flow direction were mounted repeat- 
edly on the lower and upper walls. The V-baffled channel flow was 
found to be fully developed periodic flow and heat transfer profiles 
at about x/D x 8 downstream of the tested channel inlet. It was 
reported that the P-vortex flow caused by the V-baffle induced 
impingement/attachment flows on the channel walls leading to 
drastic increase in the heat transfer rate. The heat transfer in the 
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V-baffled channel was found to be about 1-21 times higher than 
the smooth channel with no baffle. However, this also results in 
friction loss ranging from about 1.1 to 225 times above the smooth 
channel. The computational domain is shown in Fig. 47. 

A numerical investigation had been carried out by Promvonge 
et al. [60] to examine laminar flow and heat transfer character- 
istics in a three-dimensional isothermal wall square channel with 
45°-angled baffles as shown in Fig. 48. In order to generate a pair 
of mainstream wise vortex flows through the tested section, 
baffles with an attack angle of 45° were mounted in tandem and 
inline arrangement on the lower and upper walls of the channel. 
Effects of different baffle heights on heat transfer and pressure loss 
in the channel were studied and the results of the 45° inline baffle 
were also compared with those of the 90° transverse baffle and the 
45° staggered baffle. The presence of 45° baffles creates P-vortex 
flow which helps to induce impingement flows on the BLE side- 
wall and wall in the baffle cavity leading to drastic increase in heat 
transfer in the square channel. The enhancement was in the order 
of 150-850% for using both the 45° baffles with BR=0.05-0.3. On 
the other hand, the heat transfer augmentation is associated with 
enlarged pressure loss ranging from 2 to 70 times above the 
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Fig. 48. (a) Inline and (b) staggered arrangement of baffles. 
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Fig. 49. Channel geometry and computational domain with baffles. 
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Fig. 50. Double pass solar air heater with fins and baffles. 
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Fig. 51. Arrangement of Z shaped baffles in a test section. 
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Fig. 52. Channel with porous baffles. 


smooth channel. They investigated that for the 45° baffles, the 
heat transfer enhancement was around 100-200% higher than that 
for the 90° baffle whereas the friction loss can be reduced at about 
10-150%. 

A numerical investigation had been carried out by Kwankao- 
meng et al. [61] to study laminar flow and heat transfer character- 
istics in a three-dimensional isothermal wall square-channel with 
30° staggered angled-baffles. They studied the effects of different 
baffle heights at a single pitch ratio (PR=3) on heat transfer and 
pressure loss in the channel. The vortex flows was created by using 
the 30° angled baffles help to induce impingement flows on the 
sidewall and the wall in the inter baffle cavity leading to drastic 
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Fig. 53. Aluminum foam structure. 


increase in heat transfer in the channel. The order of enhancement 
is about 100-650% for using the angled baffles with BR—0.1-0.3. 
The heat transfer augmentation was associated with enlarged 
pressure loss ranging from 1 to 17 times above the smooth 
channel. Thermal enhancement factors for the inclined baffles at 
BR— 0.15 was found to be highest about 2.9. The channel geometry 
with computational domain is shown in Fig. 49. 

The performance of a solar air heater featured with double-pass 
as well as fins and baffles design was investigated by Chii-Dong 
et al. [62] using both experimental and theoretical approaches, for 
the effect of recycling operation. The performances of the new 
design and the other two designs, which are double-pass without 
fins and double-pass with fins only, were compared. The best 
collector efficiency was reported for new design as compared to all 
other designs under recycling operation. But mainly, the improve- 
ment of collector efficiency was found due to use of baffles. The 
best possible reflux ratio for the fined plus baffled double-pass 
design was investigated about 0.5 while considering both the 
collector efficiency and the pumping power requirement. A sche- 
matic drawing of a fined plus baffled double pass solar air heater is 
shown in Fig. 50. 

The influence of baffle turbulators on heat transfer augmenta- 
tion in a rectangular channel of aspect ratio of 10 fitted with the 
in-phase and out-phase 45° Z-baffles in the turbulent regime from 
Re=4400 to 20,400 had been investigated experimentally and 
numerically by Sriromreun et al. [63]. In the experiment, the 
baffles were placed in a zigzag shape (Z-shaped baffle) aligned in 
series on the isothermal-fluxed top wall, similar to the absorber 
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plate of a solar air heater channel. The intension of using the 
Z-baffles was to create co-rotating vortex flows having a signifi- 
cant influence on the flow turbulence intensity leading to higher 
heat transfer enhancement in the tested channel. Effects of the 
Z-baffle height and pitch spacing length were examined to find the 
optimum thermal performance for the Reynolds number from 
4400 to 20,400. The Z-baffles inclined to 45* relative to the main 
flow direction were characterized at three baffle to channel-height 
ratios (e/H=0.1, 0.2 and 0.3) and baffle pitch ratios (P/H— 1.5, 
2 and 3). 

For the in-phase baffle array at P/H— 1.5, the increases in Nu 
with the e/H=0.1, 0.2 and 0.3 were in the range of 430-4405, 530- 
550% and 640-670% over the smooth channel, respectively. More- 
over, the use of Z-baffles with e/H—0.3 were reported to give 
higher heat transfer than that with e/H—0.2 and 0.1 around 10% 
and 30%, respectively. The friction factor value of the Z-baffle at 
P/H=1.5 and e/H=0.3 was found to be around 2 and 5 times 
higher than that with e/H—0.2 and 0.1, respectively. 

For the in phase Z-baffle at e/H=0.1, the increases of the Nu 
for P[H— 1.5, 2 and 3 were in a range of 430-4405, 310-380% and 
310-3305: over the smooth channel, respectively. Importantly, the 
Z-baffle with e/H=0.1 at P/H— 1.5 shows higher heat transfer rate 
than the one at P/H—2 and 3 around 15% and 30%, respectively. 

The arrangement of baffles is shown in Fig. 51. 


5.19. Porous and perforated baffles 


Numerical predictions on the turbulent fluid flow and heat 
transfer characteristics for rectangular channel with porous baffles 


Balsa wood 


1/16" 


which are arranged on the bottom and top channel walls in a 
periodically staggered way was studied by Yang and Hwang [64]. 
The Experimental parameters studied include the entrance 
Reynolds number Re (1 x 104-5 x 10^), the baffle height (h/H=0.25, 
0.5 and 0.75 mm), porosities (0—0.2, 0.42 and 0.7) and kind of 
baffles (solid and porous); whereas the baffle spacing S/H are fixed 
at 1.0 and the working medium used by them was air. They found 
that if the solid baffles were replaced by porous type baffles, there 
was reduction in effective thermal conductivity of the baffle, heat 
transfer surface was increased and the flow transport phenom- 
enon was changed. Friction factor had a lower value for the 
porous-type baffle channel relative to the solid-type baffle channel 
because of less channel blockage. A schematic diagram of forced 
convection enhancement in a channel using porous baffles is 
shown in Fig. 52. 

Ko and Anand [65] carried out experimental investigation to 
measure module average heat transfer coefficients in uniformly 
heated rectangular channel with wall mounted porous baffles. 
They mounted baffles on top and bottom of the walls alternatively. 
Investigation reported that the heat transfer enhancement ratio 
decreases with increase in Reynolds number and increases with 
increase in pore density. The heat transfer enhancement ratio was 
found to be higher for taller (Bj/Dj — 2/3) and thicker (BiDg— 1/3) 
baffles. The friction factor slightly decreased with increase in 
Reynolds number, and increased with baffle thickness and pore 
density. Material used for the baffles are shown in Fig. 53. Details 
of the test section are shown in Fig. 54. 

Experimental investigations were carried out by Tzeng [66] to 
determine the local and average heat transfer characteristics in 
asymmetrically heated sintered porous channels with metallic 
baffles. The solid baffles were inserted periodically into the 
sintered metallic materials in four modes as shown in Fig. 55. 
They also examined the effect of the bead diameter of sintered 
metallic materials. The heat transfer by convection in all modes 
was increased as the bead diameter declined. It was studied that at 
Re > 2000, heat transfer enhancement was found to be more in 
mode B and worst in mode D even less than that in mode A. At a 
Re of about 1000, there was heat transfer enhancement in mode D. 


5.20. Perforated baffles 


Lin [67] analyzed the experimental results of heat transfer for 
baffles designed with different heights (H— 10-50 mm) and pores 
(N=1-3) in the event of five Reynolds numbers (702-1052) and 
three heating quantities (qi; — 90-750 W/m?) as well as heat flux 
(Q—40-100 l/min). It was studied that maximum Nusselt number 
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Fig. 55. Modes of test specimens and measured positions of thermocouples. 
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without baffles of pores was generated at qi; — 750 W/m’, Re— 1752 
at X/L— 0.04. It was found that the distribution of Nusselt number 
declines gradually with the increase of X/L, and its minimum value 
was generated at X/L —0.68, with a nearly one time difference with 
the maximum value. The biggest difference occurred at X/L— 0.2, 
with the difference of Nusselt number approx. 50%. In the case of 
XIL > 0.68, Nusselt number gradually increases, so Nusselt number 
at rear section is higher than that at middle section, i.e. outlet 
effect. 

Experimental study was performed by Huang et al. [68] for a 
measurement of heat transfer coefficients in a square channel with 
a perforation baffle by using the transient liquid crystal technique, 
and also studied the effects of Reynolds number and the height of 
the perforation baffle on Nusselt numbers. Local heat transfer 
enhancement at the downstream of step baffle was found to be 
more significant when Reynolds number was large and when 
baffle height was higher. Heat transfer coefficient off center was 
found to be better because of secondary flows that appeared off 
center after the airflow passes through the baffle. The heat transfer 
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Fig. 56. Fully perforated baffles and half perforated baffles with two and single 
rows of holes respectively. 
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Fig. 57. Different types of longitudinal vortex generators. 


vonex cae 


4 84 9^ 
A! IX ri "E 


AAA 


Flow 


Vortex Generators m 


enhancement in case of a baffle with holes was reported greater 
than that without holes. 

Karwa and Maheshwari [69] made an experimental study of 
heat transfer and friction in a rectangular section duct with fully 
perforated baffles (open area ratio of 46.8%) or half perforated 
baffles (open area ratio of 26%) at relative roughness pitch of 7.2- 
28.8 affixed to one of the broader walls as shown in Fig. 56. They 
investigated that there was an enhancement of 79-169% in Nusselt 
number over the smooth duct for the fully perforated baffles and 
133-274% for the half perforated baffles while the friction factor 
for the fully perforated baffles was found to be 2.98-8.02 times 
of that for the smooth duct and 4.42-17.5 times for the half 
perforated baffles. 


5.21. Delta winglet 


In order to enhance the rate of heat transfer in engineering 
devices like heat exchanger, vortex combustor, drying process etc, 
secondary or rotary flow is created using swirl or vortex gen- 
erators that are inserted in main flow. These are made in triangular 
or rectangular shapes and can be welded or punched and bend out 
of the plate so that they project into the flow with an angle of 
attack to the main flow direction as shown in Fig. 57. Using delta 
wing vortex generator, enhancement in average heat transfer at 
low Reynolds number was found to be about 50%-60% [70]. 

An experimental study was carried out by Torii et al. [71] to 
obtain heat transfer and pressure loss in a test section, simulated 
to a fin-and-tube heat exchanger, with in-line or staggered tube 
banks with delta winglet vortex generators of various configura- 
tions as shown in Fig. 58. The proposed "common flow up" config- 
uration causes significant separation delay, reduces form drag, and 
removes the zone of poor heat transfer from the near-wake of the 
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Fig. 59. Schematic of delta winglet. 
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Fig. 58. Configuration of winglet type vortex generator on the fin surface-tube bank: (a) “common flow down" configuration; (b) “common flow up" configuration. 
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tubes. In case of staggered tube banks, the heat transfer was 
augmented by 30% to 10%, and yet the pressure loss was reduced 
by 55% to 34% for the Reynolds number (based on two times 
channel height) ranging from 350 to 2100, when the present 
winglets were added. In case of in-line tube banks, these were 
found to be 20% to 10% augmentation, and 15% to 8% reduction, 
respectively. 

An experimental study was carried out by Yakut et al. [72] to 
study the effects of various kinds of design parameters on the heat 
transfer and flow characteristics of the tapes with delta-winglet 


a b 


Top View 


vortex generators were systematically analyzed using the Taguchi 
method. The maximum heat-transfer rate was observed at 16,906 
Reynolds number, 60° angle of attack, 16 mm winglet height and 
25 mm pitch of the winglets and the minimum friction factor was 
observed at 8 mm winglet height, 16906 Reynolds number, 30* 
angle of attack and 75 mm pitch. The maximum vortex-shedding 
frequency occurred at 16,906 Reynolds number, 25 mm pitch, 
8mm winglet height and 30* angle of attack. The minimum 
vortex-shedding frequency was obtained at 3690 Reynolds num- 
ber, 75 mm pitch, 16 mm winglet height and 90* angle of attack. 
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Fig. 60. (a) Winglet vortex generator dimensions; (b) single row; (c) three row inline array configurations. 


Staggered arrangement 


In-line arrangement 


Fig. 61. Schematic of delta winglets in (a) staggered and (b) inline arrangement. 


Fig. 62. Schematic view of present modified longitudinal vortex generator: (a) layout of common-flow-down wings; and (b) modified rectangular wing. 
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The optimum results were obtained at 25 mm pitch, 8 mm height 
of winglet, 30° angle of attack and a Reynolds number of 16,906 
for larger frequencies required to avoid resonance. For smaller 
frequencies, in order to avoid resonance, the optimum results were 
75 mm pitch, 8 mm height of winglet, 60* angle of attack and a 
Reynolds number of 16,906. Geometry used for the investigation is 
shown in Fig. 59. 

The potential of winglet type vortex generator (VG) arrays for 
air-side heat transfer enhancement was experimentally evaluated 
by full-scale wind-tunnel testing of a compact plain-fin-and-tube 
heat exchanger by Joardar et al. [73] as shown in Fig. 60. They 
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studied that the air-side heat transfer coefficient increases from 
16.5% to 44% for the single-row winglet arrangement with an 
increase in pressure drop of less than 12%. For the three-row 
vortex generator array, the enhancement in heat transfer coeffi- 
cient increases with Reynolds number from 29.9% to 68.8% with a 
pressure drop penalty from 26% at Re=960 to 87.5% at Re=220. 

Tian et al. [74] carried out a three-dimensional numerical 
simulation to study the air-side performance of a wavy fin-and- 
tube heat exchanger with delta winglets in terms of heat transfer 
and fluid flow characteristics. The wavy fin-and-tube heat exchan- 
gers which have three-row round tubes in staggered or inline 
arrangements were studied. They studied that when , Rep-=3000, 
compared with the wavy fin, the j and f factors of the wavy fin 
with delta winglets in staggered and in-line arrays are increased 
by 13.1%, 7.0% and 15.4%, 10.5%, respectively. Geometry investi- 
gated is shown in Fig. 61. 

Fluid flow and heat transfer characteristics of a modified 
rectangular longitudinal vortex generator (LVG) as shown in 
Fig. 62, obtained by cutting off the four corners of a rectangular 
wing was studied by Min et al. [75] experimentally. Study shows 
that the modified rectangular wing pairs (MRWPs) have better 
flow and heat transfer characteristics than those of rectangular 
wing pair (RWP). It was found that at z= + 40 mm from the 
centerline of the heater plate, the local heat transfer is enhanced 
due to the strong longitudinal vortices generated by the presence 
of the LVGs. 

Second law analysis of a cross-flow heat exchanger (HX) 
was done by Kotcioglu et al. [76]. The entropy generation in a 
cross-flow HX with a new winglet-type convergent-divergent 


Top view 


Fig. 64. Twisted tape vortex generator: (a) typical twisted tape (TT), (b) straight delta-winglet twisted tapes (S-DWT) and (c) oblique delta-winglet twisted tapes (O-DWT). 
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longitudinal vortex generator (CDLVG) is investigated and is 
shown in Fig. 63. They found that at the low Reynolds number, 
the entropy generation number was influenced by the heat 
transfer while at the high Reynolds number, it was influenced by 
the pressure drop. Another conclusion reported that, the present 
vortex generator (CDLVG) shows an increase in the heat transfer 
enhancement from 15% to 30% and also an increase in the 
pressure-loss penalty from 2055 to 3055, in a comparison with and 
without vortex generators, respectively. 

Experimental investigations of heat transfer, flow friction and 
thermal performance factor characteristics in a tube fitted with 
delta winglet twisted tape, using water as working fluid was done 
by Eiamsa-ard et al. [77]. In the study, they described the 
influences of the oblique delta-winglet twisted tape (O-DWT) 
and straight delta-winglet twisted tape (S-DWT) arrangements. 
They reported that the O-DWT is more effective turbulator giving 
higher heat transfer coefficient than the S-DWT. Over the range 


considered, Nusselt number, friction factor and thermal perfor- 
mance factor in a tube with the O-DWT are, respectively, 
1.04-1.64, 1.09-1.95, and 1.05-1.13 times of those in the tube with 
typical twisted tape (TT). Geometry investigated is shown in 
Fig. 64. 

Champookham et al. [78] carried out experimental investiga- 
tions to study the effect of combined wedge ribs and winglet type 
vortex generators (WVGs) on heat transfer and friction loss 
behaviors for turbulent airflow through a constant heat flux 
channel. Two types of wedge rib viz. wedge ribs pointing down- 
stream and pointing upstream were arranged in the form of inline 
and staggered style. Two pairs of the WVGs with the attack angle 
of 60° were mounted on the test channel entrance so as to 
generate longitudinal vortex flows through the tested section. 
The test channel had an aspect ratio, AR=10 and height, 
H=30mm with a rib height, e/H—0.2 and rib pitch, P/H— 1.33 
with the Reynolds number range of 5000-22,000. Investigations 
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Fig. 65. (a) winglet geometry, (b) wedge rib pointing upstream and (c) wedge rib pointing downstream. 
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Fig. 66. Delta winglets pointing downstream (PD) and pointing upstream arrangement. 
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Fig. 67. Different types of vortex generators under study. 
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reported that the use of the wedge rib turbulators with e/H=0.2 
causes a high pressure drop increase, specially for the in-line 
wedge pointing downstream and also provides considerable heat 
transfer augmentations, Nu/Nuo— 2.9-3.5 for use with the WVGs. 

The geometry investigated is shown in Fig. 65. 

Promvonge et al. [79] experimentally studied the effects of 
combined ribs and delta-winglet type vortex generators (DWs) on 
forced convection heat transfer and friction loss behaviors for 
turbulent airflow through a solar air heater channel. Rectangular 
channel with aspect ratio, AR=10, height, H—30 mm with 
Reynolds number in the range of 5000-22,000 was used. They 
reported that combined rib and the DW provides considerable 
heat transfer augmentations, Nu/Nuo—2.3-2.6 and also causes a 
moderate pressure drop increase, f/fo— 4.7—10.1, depending on the 
attack angle and Re values. Another conclusion is that, the use of 
combined rib and PD-DW at lower angle of attack provides higher 
heat transfer of about 40-65% and better thermal performance 
than the rib/the DW alone, leading to more compact heat exchan- 
ger. Delta winglet geometries are shown in Fig. 66. 

The performance of a pair of new vortex generators called 
Curved trapezoidal winglet (CTW) had been investigated experi- 
mentally by Zhou and Ye [80] and compared with traditional 
vortex generators - rectangular winglet, trapezoidal winglet and 
delta winglet as shown in Fig. 67. They studied that Delta Winglet 
Pair (DWP) is the best in laminar and transitional flow region, 
while Curved Trapezoidal Winglet Pair (CTWP) has the best 
thermo hydraulic performance in fully turbulent region due to 
the streamlined configuration and then the low pressure drop, 
which indicates the advantages of using this kind of vortex 
generators for heat transfer enhancement. Small attack angle of 
CTWP, such as fj—-0^ and 15°, have better thermo hydraulic 
performance of R than larger attack angles. Particularly, when Re 
is larger than 18,000, /—0* presents the best thermo hydraulic 
performance with R as high as 1.6. 


5.22. Use of impinging jet 


Experimental investigation has been carried out by Chauhan et 
al. [85] to study heat transfer and friction factor characteristics 
using impinging jets in solar air heater duct. The study shows that 
there is considerable enhancement in heat transfer and friction 
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factor by 2.67 and 3.5 times respectively. Impingement plate 
geometry is as shown in Fig. 68. 


6. Computational analysis 


Literature survey in the field of artificially roughened solar air 
heater reveals that, there are very few attempts had been made in 
early days because of hardware and development in numerical 
methodology to solve the complex flow equations. A CFD based 
investigation of turbulent flow through a solar air heater is carried 
out by Yadav et al. [86] with square sectioned transverse rib 
roughness as shown in Fig. 69. The 2D analysis of rectangular duct 
is performed using ANSYS FIUENT 12.1 software and using RNG k-e 
turbulence model. Analysis using four different configurations of 
rib roughness keeping relative roughness pitch constant at P/e 
— 14.29 and six different values of Reynolds number, ranging from 
3800 to 18,000, reveals that the relative roughness height is a vital 
factor and mainly affects the rate of heat transfer and the increase 
in flow friction. In this analysis, the maximum value of thermo 
hydraulic performance parameter (THPP) has been found to be 
1.8 for the range of parameters investigated. 

Another numerical investigation is conducted by Yadav et al. 
[87] using FLUENT 12.1 and RNG k-e turbulence model on the 
solar air heater, using square sectioned transverse ribs as shown in 
Fig. 70, with different values of relative roughness pitch (7.14 « P/e 
x 17.86). Reynolds numbers are in the range of 3800 to 18000. 
This study reported the value of THPP as 1.82 by the use of ribs 
with P/e of 10.71. 

2D CFD analysis to predict the effect of small diameter of 
transverse wire ribs roughness on heat transfer enhancement in 
solar air heater is done by Yadav et al. [88]. In this analysis also 
FLUENT 12.1 software is used with the RNG k-e turbulence model. 
Here the operating parameters are taken as, (P/e=7.14-35.71) and 
(e/D=0.021-0.042). Here the maximum value of THPP is found to 
be 1.65 for the range of parameters investigated. The domain used 
for CFD analysis is shown in Fig. 71. 

CFD investigation using ten different ribs namely rectangular, 
square, chamfered, triangular etc. is done by Chaube et al. [89]. 
FLUENT CFD code is used for analysis using the SST k-e turbulence 
model. The heat flux of 1100 W/m? is provided on the absorber 
plate only and rib surface is kept adiabatic. Their study reported 
that highest heat transfer is achieved with chamfered ribs and the 
best performance index is found with rectangular ribs of size 3 x 5. 
Two dimensional domain is shown in Fig. 72. 

A three dimensional (3D) solar air heater duct provided with 
artificial roughness in the form of thin circular wire in arc shaped 
as shown in Fig. 73, has been analyzed by kumar et al. [90] using 
CFD. Renormalization (RNG) k-e turbulence model is used for 
analysis in FLUENT CFD code. Their study reveals the maximum 
value of overall enhancement ratio as 1.7 for the range of 
parameters investigated. 


Square transverse ribs 


> 


ir out 


i 


= 20mm—>! 


k—H 


Fig. 69. Geometry of two dimensional domain. 
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Fig. 71. Two dimensional domain for CFD analysis. 
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Fig. 72. (a) Ribbed duct geometry for CFD analysis and (b) two dimensional computational domain. 


Fig. 73. Geometry of roughened surface collector plate. 


Karmare et al. [91] had carried out analysis using CFD on solar 
air heater with a absorber plate roughened with metal ribs of 
circular, square and triangular cross-section (Figs. 35 and 74), 
having 54^, 56°, 58°, 60° and 62° inclinations to the air flow. 
The system and operating parameters studied were: e/D— 0.044, 
P/e=17.5 and l/sc 1.72. Range of Reynolds number was 3600- 
17000. Their study reported that maximum heat transfer is 
obtained with the square cross-section ribs with 58 angle of 
attack and there is 30% enhancement in the heat transfer for 
square plate over smooth surface. Fig. 74. Flow domain with 60° rib having circular cross section. 
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Fig. 75. Comparison of enhancement factor of different geometries. 


The effect of wedge-shaped transverse ribs roughness on flow 
through a rectangular duct of a solar air heater is studied using 
CFD by Gandhi et al. [92]. Two dimensional CFD simulation is 
carried out using FLUENT on a duct having relative roughness 
pitch (P/e) of 4.5, relative roughness pitch of 0.022 and rib wedge 
angle of 15*. The results of numerical analysis are in good 
agreement with the experimental results in terms of velocity 
profile and turbulence intensity. The other observations revealed 
that rib roughened surface in a rectangular duct produced much 
flatter velocity profile and higher turbulence intensities as com- 
pared to smooth surface. This causes increase in heat transfer in 
roughened duct. Friction factor for roughened duct was around 
3.4 times higher than the smooth duct at Reynolds number of 
21,000. 

CFD analysis on solar air heater is done by Sharma et al. [93] using 
V-shaped ribs roughness on the underside of the absorber plate. 
Finite volume method with semi-implicit method for pressure linked 
equations (SIMPLE) algorithm is used for computations. Reynolds 
number range was from 5000 to 15,000. The results obtained shows 
that combined effect of swirling motion, detachment and reattach- 
ment of the fluid were responsible for the increase of heat transfer 
rate during CFD analysis. Nusselt number increases and friction factor 
decreases with increase in Reynolds number for all combination of 
relative roughness height (e/D) and relative roughness pitch (P/e). An 
average percentage deviation predicted between CFD and exact 
solution was found less than - 34. V-shaped rib roughness found 
to give high rate of heat transfer. 


7. Comparison of thermo hydraulic performance of roughened 
solar air heaters 


Various researchers reported that there is a heat transfer 
enhancement because of the inclusion of the roughness geometry 
in the path of the air flow. But this heat transfer enhancement is 
also leads to increased pumping power penalty due to correspond- 
ing increase in friction factor. So it is essential to determine the 
geometry that will results in the maximum enhancement in heat 
transfer and minimum increase in friction factor. This is done 
using a factor developed by Lewis [94] known as efficiency 
parameter, 7 which evaluates the enhancement of heat transfer 
for same pumping power requirement. It is defined as 


q= (Nu, /Nuj)/(, /f ) ? 
Correlations of heat transfer and friction factor for different 
roughness geometries are used to derive thermo hydraulic 


performance. Comparisons are made for the geometries having 
relative roughness pitch (P/e) value as 10 and relative roughness 
height 

(e[D) value as 0.04. Range of Reynolds number is taken in between 
2000 and 20,000. Thermo hydraulic performance parameter 
comparision of different types of ribs has been shown in Fig. 75. 
Following conclusions may be drawn from the figure 


(1) Amongst all the roughness geometries considered for the 
comparison, Multi V-rib with a gap is found to give best 
thermohydraulic performance. 

(2) The next best performance after multi V-rib with gap is 
achieved by V-rib as shown in Fig. 75. 

(3) Other ribs are found to give thermo hydraulic performance in 
the range of 0.5 to 2.0. 


8. Conclusions 


In the present paper, an attempt has been made to review 
roughness geometries which have been used by researchers for 
the enhancement of heat transfer in solar air heater. It is found 
that substantial enhancement in the heat transfer can be achieved 
with little penalty of friction. In general the salient features of the 
reviews are 


(1) Many investigators carried out number of experimental and 
numerical studies for better understanding of heat transfer 
enhancement and flow processes due to turbulence generated 
by the presence of roughness on absorber plate of a solar air 
heater. The paper presents a review of experimental as well as 
numerical analysis carried out by the researchers. 

(2) The relative roughness height, relative roughness pitch, angle 
of attack and relative gap position are some of the parameters 
studied for their effect on heat transfer and friction character- 
istics in both heat exchangers and solar air heaters. The paper 
has given geometrical parameters for individual geometry, 
specifying the heat transfer enhancement in terms of Nusselt 
number and pressure drop penalty in terms of friction factor 
values. 

(3) The application of solar air hater varies from industry to 
industry. Depending on the energy requirement, particular 
type of geometry can be selected. The reviews given in this 
paper may be useful in this regards. 

(4) The heat transfer enhancement and friction factor character- 
istics in solar air heater are studied recently using Computa- 
tional Fluid Dynamics (CFD) i.e. ANSYS-FLUENT software. CFD 
simulations results are found to be in good agreement with 
experimental results. The literature review on this area pre- 
sented in this paper along with experimental work will be 
beneficial for the researchers working in simulation area 
using CFD. 

(5) The correlations based on experimental results in the form of 
Nusselt number and friction factor are available in earlier 
papers which may be used to predict the thermal as well as 
thermo hydraulic performance of solar heaters. Heat transfer 
coefficient and friction factor correlations reported in the 
literature are presented in Table 1 which may be useful to 
the beginners for doing research in this area. 

(6) Comparisons between various roughness geometries can be 
made on the basis of a parameter known as thermo hydraulic 
performance parameter (THPP). It is observed that multi V-rib 
with a gap gives better performance in comparison to others 
roughness geometries for Reynolds number range between 
2000 and 20,000. 


Table 1 
Heat transfer coefficient and friction factor correlations for different roughness geometries used in solar air heater duct. 


Investigators Year Roughness geometry Range of parameters Correlations 
Heat transfer Friction factor 

Prasad and Saini [8] 1988 Transverse rib T ae S = (f/2)/1+ Vf/245(e* y?28pr957 _ Q.95(P/e)°53y] f =2/[0.95(P/e)°53 42.5 In(D/26) 3.75] 

Re: 5000-50,000 
Lau et al [27] 1991 Combined rib (Transverse, = G(e* ) = a(e* a = 1.768 — 3.884b = 0.229 — 0.333 R(e*) = a(e* )"a = 0.968 — 2.269b = 0.055 — 0.152 

inclined and V-shaped rib) 

Hong and Hsieh [82] 1991 Transverse ribs e/Dp: 0.19 Nu-0.056Re,97^ (for Tu — 0.04) = 

ple: 5.31 Nu-0.030Re497? (for Tu = 0.07) 


Nu=0.019Re,°? (forTu- 0.11) 


Gupta et al.[18] 1994 Transverse ribs e[D: 0.018-0.052 Nu = 0.000824(e/D) - 78 (Ww /Hy?84Re!062? e —35 f — 0.06412 (e/D)9 99 w /H)0237 Re - 0185 
Re: 3000-18,000 Nu = 0.00307(e/D) "^9 (w/Hy Re. e 35 

Gupta et al. [22] 1997 Inclined Ribs e[D: 0.02-0.053 Nu = 0.000824(e/D) - 178 (W /H 9284 Re!.062 f —0.06412(e/D)9919 (W /Hy0237 Re - 0.185 
ple: 7.5-10 exp[—0.04(1—a/60)?(k/D)e < 35 exp[—0.0993(1 —2/60)?] 
a : 30-90* 


Nu= 0.00307(e/D)~°489(W/H)2245 Re0812 
Re: 5000-30,000 
exp[—0.475(1—a/60)^k/D)e > 35 
Saini and Saini [40] 1997 Expanded wire mesh metal ribs lle: 25-71.87 Nu = 4.0 x 107 ^Re! 2(e/Dj995(s 19g??? f = 0.815Re9?6! (/g)0266(5 /1 9g)919 (10e /d)?59! 
sle: 15.62-46.87 exp[- 1.25( In(s/106))?] x (1/10e)59 


e[D: 0.012-0.0390 2 
Re: 1900-13,000 xexp[—0.824{ In (1/10e)*}] 


Muluwork et al.[33,34] 1998 V shaped staggered e[D: 0.02 Nu = 0.00534(B/S) ?^95ge! 2991 f 5 0.7117 (B/S 9636gg 2991 
Discrete wire ribs a: 60* 
B/S: 3-9 
Re: 2000-15,500 
Karwa et al. [39] 1999 Chamfered ribs e[D: 0.014-0.0320 G = 103.77e - 9906 (w /H)95 (p /e^56 for 7«xe* «20 
ple: 4.5-8.5 exp[0.7343( In(p/e))?e* ) ?! for 7 xe* «20 R— 1.66e - 00078 (uw /]) -04 (p /gy2.695 
Re: 3000-20,000 G — 32.26e 995 (w /Hy95(p/ey 56 exp[—0.762{ In(p/e)}2(e+ )- 9°75] 
@: — 15, 0, 10, 15, 18 2 08 + 
exp[0.7343{ In(p/e))?y(e * )°* for 20«e* <60 for 20<e* <60 


WIH: 4.8, 6.1, 7.8, 9.66, 12 
R= 1.325e - 00078 (w /H)~°4(p/e)?695 


exp[— 0.762( In(p/e)}?] 


Varma and Prasad [19] 2000 Transverse protrusion wire e[D: 0.01-0.03 Nu = 0.08596(p/e) - 995^ (e /D)9.972 ge0723e < 24 f —0.0245(p/e) 9 9?06(e 0.921 gg - 125 
ple: 10-40 Nu = 0.0245(p/e) - 916 (e /p)90?1 Re? 802¢ > 24 
Re: 5000-20,000 
e*:8-42 
P,: 0.7 

Momin et al.[31] 2002 V-shaped ribs e[D: 0.02-0.034 Nu = 0.067(e/D)9 4^ (a /60)7 9.077 Re 888 f — 6.266(e/D)9595 (a/60) - 0993 Re - 0425 
ple: 10 exp[— 0.782 In (a/60)?] exp[—0.719 In (a/60)?] 
a: 30-90* 


Re: 2500-18,000 


Bhagoria et al [36] 2002 Wedge shaped ribs e[D: 0.015-0.033 Nu = 1.89 x 107 ^Re!?! (e/D)94?6 (p/e) (p / 9) - 0.018 f —1244(e/DY 99 (6/10)? 49Re -.18(5 /g) -052 


Bs [exp(— 0.71 n(p/e))?)Jfexp(— 1.50n(9/10)?)] 
a. lcd! D 


Re: 3000-18,000 


PLE 
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Karwa [35] 


Jaurkeret al. [37] 


Karmare and Tikewar 
[46] 


Layek et al.[38] 


Saini and Verma [43] 


Saini and Saini [45] 


Kumar et al. [48] 


Varun et al.[26] 


Bopche and Tandale [49] 


Aharwal et al. [23] 


Karwa et al. [69] 


Hans et al. [11] 


2003 


2006 


2007 


2007 


2008 


2008 


2008 


2008 


2009 


2009 


2009 


2010 


Inclined discrete and continuous wire 
ribs 


Rib-groove 


Metal grit rib 


Chamfered rib groove combination 


Dimple shaped protrusion 


Arc shaped wire ribs 


Discrete W shaped ribs 


Inclined and transverse wire ribs 


Inverted U shaped solid baffles 


Inclined discrete ribs 


Transverse perforated baffles 


Multiple V shaped 


e[D: 0.0467-0.05 
ple: 10 

a: 60-90* 

BIS: 3 

WIH: 719-7.75 


e[D: 0.0181-0.0363 
ple: 4.5-10 

glp: 0.3-0.7 

Re: 3000-21,000 


e[D: 0.035-0.044 
ple: 12.5-36 

l[s: 1.72-1 

Re: 4000-17,000 


e[D: 0.022-0.04 
ple: 4.5-10 

glp: 0.3-0.6 

Re: 3000- 
21,000 

e[D: 0.018-0.037 
ple: 8-12 

Re: 2000-12,000 


e/D: 0.0213-0.0422 
ple: 10 

gle: 0.5-2 

a/90: 0.333-0.666 
glp: 0.3-0.6 

Re: 2000-17,000 


e[D: 0.0338 
ple: 10 ° 
a: 60 


Re: 2000-14,000 
ple: 3-8 


e/D: 0.0186-0.03986 
ple: 6.67-57.14 

a: 90° 

Re: 3800-18,000 


e/D: 0.0377 

ple: 10 

WJH: 5.87 

d/w: 0.167-0.5 
a: 60° 

Re: 3000-18,000 


Re: 2700-11,150 

ple: 7.21, 14.42, 28.84 
B: 26%, 46.8% 

d/e: 0.047 


e/D: 0.019-0.043 
ple: 6-12 

a: 30-75* 

Re: 2000-20,000 
W/w: 1-10 


Cs 32.26e - 9906 (W /H)95 (p /e)?56 
exp[0.7343{ In(p/e)? (e * ) 998 


Nu = 0.00206Re^935 (e /)9349(p 3318 
xexp[—0.868{ In(p/e))?g/p). 98 
xexp[2.486( In(g/p))? + 1.405{ In(g/p))?] 


Nu= 2.4Re!/3(e/D)°*?(I/s)~°-48(p/e) -027 


Nu = 0.002062Re 939 (e / 0349 (p /gj3318 
exp[—0.868{ In(p/e)? (g/p) 18 
xexp[2.486( In(g/p))? + 1.406( In(g/p))?] 


Nu =5.2 x 10> ^Re! 7 (e/D)93? (p/ey-15 


exp[— 2.21( In(p/e)?] x exp[- 1.30( In(e/D))?] 


Nu = 0.001047(e/Dy?3772(a/90) - 0:198 ge 13186 


Nu— 0.105Re^57 (p; 9453 ( /g0) - 0 081 
exp[—0.59 x (In(a/60))?] 


Nu = 0.0006(p/e)9 9104 ge! 213 


Nu= 0.5429Re? 7054 (@ / [) 0.3619 (p jg) - 0 1592 


Nu = 0.0102Re! !45(e /pj9^! 
[(1—(0.25 — d/wy (0.0101 —g/e)?)}] 


Nu = 0.0893Re? 7608 


Nu —3.35 x 10 5Re?9?(e/D9 "7 (W /w)93 (q/90) - 949 
exp[—0.1177(In(W /w))?Jexp[ — 0.61 (In(a/90))?] 


(p/e^^exp[— 2.0407(n(p/e)?] 


for 7xe* «20 
R= 1.66e- 9-078 (W=/H)~°4(p/e)?995 


exp[—0.762{ In(p/e)}?(e* ) 9975] 
for 20<et «60 

R= 1.325e 9-078 (w /H) -94(5 Jg? 695 
exp[ — 0.762( In(p/e)}?] 


f = 0.0012Re 9:199 (e /)9585 (p /g) 19 
xexp[— 1.854( In(p/e))?(g/p)? ^^ 
xexp[1.513( In(g/p))? +0.8662{ In(g/p)?] 


f= 15.55Re  9?9(e/ [99^ (l/s) -027 (p /g) - 051 


f 5 0.001227Re 9199 (e /)9585 (y /g)7.19 
exp[—1.854{ In(p/e))?(g/p)9 ^ 


exp[1.513( In(g/p)}? --0.8662( In(g/p)?] 


f — 0.0642Re - 942 (e) - 00214 (p jg) -0465 
exp[0.054{ In(p/e))?Jexp[ — 0.840( In(e/D))?] 


f- 0.14408(e/D)? 795 (a /90)9 1185gg — 0.17103 


f =5.68 x Re~°*°(e/D)°5°(a/60)~ 9-8! 
exp[—0.579 x (In(a/60))?] 


f = 1.0858(p/e)9! I4 Re - 03685 


f= 1.2134Re~ 02076 (e /[))0.3285 (p /e)— 94259 


f =0.5Re =00836 (@ (DOr? 


f =0.1673 Re 99215 


f =447 x 10 ^Re 03188 (e /D)°-73 (W/w)??? 
(a/90) ° exp[ — 0.52(1n(«/90))^(p/e)?9 
exp[ - 2.133 (n(p/e)?] 
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Table 1 (continued ) 


Investigators 


Promvonge [55] 


Kumar et al. [25] 


Singh et al. [32] 


Bhushan and Singh [44] 


Lanjewar et al. [47] 


Yadav et al. [83] 


Sethi et al. [84] 


Chauhan et al [85] 


Year 


2010 


2011 


2011 


2011 


2011 


2012 


2012 


2013 


Roughness geometry 


Multiple 60* V shaped baffles 


60° inclined continuous discrete rib 


Discrete V-down ribs 


Dimple Shape Roughness 


W-Shaped ribs 


Circular protrusions 


Angular fashion dimples 


Impinging jets 


Range of parameters 


e/H:0.10, 0.20, 0.30 
PJH: 1, 2,3 
Re: 5000-25,000 


Re: 4105.2-20,526.2 

e/D: 0.0249, 0.0374 and 0.0498 

p/e: 8, 12 and 16 d/W: 0.15, 0.25 and 0.35 
gle: 1 


Re: 3000-15,000 
ple: 4-12 

a: 30-75* 

d/w: 0.2-0.8 
gle: 0.5-2.0 

e[D: 0.015-0.043 


Sle: 18.75-37.50 
L/e: 25-375 
d/D: 0.147-0.367 
e/D: 0.03 

Re: 4000-20,000 


ple: 10 

e: 15 

e[D: 0.03375 
WIH: 8 

Re: 2300-14,000 
a: 30, 45, 60-75* 
e*: 8-44 


Re: 3600-18,100 
ple: 12-24 

e[D: 0.015-0.03 
a: 45-75 

WJH: 11 

eld: 0.3 


WJH: 11 

e[D: 0.021-0.036 
ple: 10-20 

a: 45-75° 

e/d: 0.5 

Re: 3600-18,000 


Re: 3800-16,000 
x|Dp: 0.435-1.739 
y[Ds: 0.435-0.869 
W/Z: 11.6 

D;|Dp: 0.045-0.109 


Correlations 


Heat transfer 


Nu = 0.147Re? 7? (pR4(1 —e/H)~ 1793 (1. pg) 94? 


Nu=3 x 10^ Re? (e 0290 (p/e)? 885 (q /w9-115 
exp[- 1.237( In(p/e))?] 


Nu =2.36 x 10 >Re? (p /g)350 (4 /gg) - 0.023 (d/w)~ 0.043 
(g/e) 99 exp[ —0.84{ In(p/e)?Jexp[ — 0.72( In(a/60))?] 
exp[—0.05{ In(d/w))?]exp[ — 0.15( In(g/e))?] 


Nu=2.1 x 107 35Re! 492 (5 /e)!2 94([/gy99 ?(d/D) -?9 
exp[— 10.4( log (S/e))?]exp[ —77.2( log (L/e))?] 
exp[- 7.83( In(d/D))?] 


R= \/(2/f)+2.5 In(2e/D) 43.75 
et =4/(f/2)Re(e/D) 


S=([(f/2St)—1}\/(2/f/)+R 


Nu= 0.154Re!°!7(p/e)~938(e/D)°*"(a/60) - 0:213 
exp[—2.023{ In(a/60))?] 


Nu= 0.154Re! 9 (p/e) -038 (e 0521 (a/60)~ 0.213 
exp[— 2.023( In(a/60))?] 


Nu — 1.658 x 10 3Re? 8572 X / Dy 0.1761 y jp, 0141 
(Dj/Dj) - 959 *exp—0.3498( In(D;/D;))?] 


Friction factor 


f = 0.48Re®-38 (1 —e/H) ??8(1 + PR)- 9833 


f = 0.0014Re~°3(e/D)° 804 (p/e)??? (p /gy*516 
exp[— 0.944( In(p/e)}?] 


f =44.13 x 10 ?Re 9 75(pjgj74(q/gQ) 994 
(d/w)~ 9958 (g jg) - 0.031 (e/D)979 

exp[—0.685{ In(p/e))?Jexp[ — 0.93( In(a/60))?] 
exp[—0.058{ In(d/w))?]exp[ — 0.21( In(g/e))?] 


f 2:32 Re - 9201 (S/e) 9383 (L/e)— 0.484 (q py0-133 


f =7.207 Re °58(p/e)~° 18(e/D)° 16 (5 /gQ)9.058 
exp[1.412{ In(@/60)}7] 


f —4.869 x 10^ 1 Re 9223 (pe) 79059 (e / [92663 Ca /60)° 0042 
exp[—0.4801{ In(@/60)}7] 


f —0.3475 Re~ °544(X/D), y 4169 (y Di 05321 
(D;/D,)~ 488 exp[ —0.2210{ In(D;/Dn)}71 


OLE 
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